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§ Hl Abstract 

The rise with energy of the neutrino-nucleon cross section implies that at energies above 
few TeV the Earth is becoming opaque to cosmic neutrinos. The neutrinos interact with the 

5-H ' 



nucleons through the weak charged current, resulting into absorption, and the weak neutral 
current, which provides a redistribution of the neutrino energy. We Mellin transform the 
neutrino transport equation and find its exact solution in the moment space. A simple 
analytical formula is provided, which describes accurately the neutrino spectrum, after the 
neutrinos have traversed the Earth. The effect of the weak neutral current is most prominent 
for an initial flat neutrino spectrum and we find that at low energies (around 1 TeV) the 
neutrino intensity is even enhanced. 
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Neutrino telescopes will provide a new observational window on our cosmos and will help 
us probe the deepest reaches of distant astrophysical objects The enormous energies 
involved (cosmic neutrinos carry energy from few TeV to 10 8 TeV) might also allow us 
to enlarge our knowledge of particle physics. Since neutrinos interact only weakly, they 
travel unhindered by intervening matter and can give us information about regions of the 
universe, which are not accessible to traditional photon astronomy. By the same token, it is 
extremely difficult to detect neutrinos and large effective volumes for detection are required. 
Four neutrino observatories are under construction: DUMAND II in the Pacific Ocean off 
the coast of Hawaii §, NESTOR in the bay of Pylos, Greece 0, AMANDA in the ice 
of Antarctica ||] and the Baikal NT-200 in the Siberian lake Baikal ||. Expected neutrino 
sources include atmospheric neutrinos (neutrinos generated in the atmosphere by the cosmic 
rays), which dominate other sources at energies below few TeV, neutrinos generated in active 
galactic nuclei (AGN) with energies up to 10 6 TeV and cosmological neutrinos, important 
at energies above 10 6 TeV 

The neutrino-nucleon cross section rises linearly with energy up to 10 TeV. At higher 
energies, due to the propagator effect of the intermediate gauge bosons, the rise is slower. At 
around 40 TeV the cross section is large enough so that the Earth starts becoming opaque to 
neutrinos. To reduce the background from cosmic rays, at the detection site we are looking 
for upward moving muons induced by neutrinos coming from the other side of the Earth. 
Therefore, in all relevant calculations we should include attenuation factors describing how 
the neutrino fluxes change as the neutrinos travel through the Earth. It is the purpose of 
the present work to provide an analytic and accurate representation of the shadowing of 
highly energetic neutrinos. 

Neutrinos interact with the nucleons through charged and neutral weak currents. In the 
first case a neutrino is transformed into a charged lepton and we have a neutrino loss, while 
in the second case the neutrino continues along its path with reduced energy. The inclusive 
neutrino-nucleon cross section looks generically like 
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d 2 a _ A(x } Q 2 ) + (l-y) 2 B(x,Q 2 ) 



dxdy " (1 + s B xyy 

In the above expression x is the Bjorken scaling variable, 1 — y is the momentum fraction of 
the produced lepton, A and S are structure functions expressed in terms of the appropriate 
parton densities and s B = s/M| with Mb the mass of the intermediate gauge boson (W or 
Z). At high energies (s B ^> 1) the propagator term (the term in the denominator) reduces 
the rise with energy of the total cross section. Also, at high neutrino energies the nucleon is 
probed at exceedingly small x values, where no experimental information is available. The 
neutrino-nucleon cross section at very high energies is a subject of current research |6|-§. 
For our purposes it is sufficient to parametrize the total cross sections in the following form 

a cc (E)=a c ( — ) , (2) 
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E \ l3n 

a nc (E) = a n ( — . (3) 



Recent calculations || provide for neutrino energies above 10 3 TeV, (3 C ~ (3 n ~ (3 ~ 0.4, 
a c /a n ~ 2.53. An important ingredient in the shadowing process is the y distribution of the 
neutral current cross section. We assume a factorized expression 

mfm. ( 4) 



dy \E , 

There is a mild energy dependence in the y distribution, but for our practical work we use 
an f(y) independent of energy. In the domain s B ^> 1 the gauge boson propagator creates 
a peaking of the cross section near y = and we adopt the form || 

f(y) = r^rr (5) 

ln(l/e) e + y 

The average y and the parameter e are related through < y >= 1/ ln(l/e). 
The transport equation for neutrinos traversing the Earth is 

= + a nc ) I(E, r) + d -^®I, (6) 



where dr = n(z)dz and n(z) is the number density of the nucleons encountered by the 
neutrino along its path through the Earth. I(E, r) is the intensity of the neutrino flux at 
"depth" r, with the initial intensity (before entering the Earth) 1(E) = I(E,t = 0). The 
convolution is defined by 



da r , 



da nc (E',y) 



/(£', t) 5(E - E'(l - y)) dE'dy. 
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dy J dy 

For n(z) we use n = p/M^, p being the density of the Earth. The maximum value r can 
attain is approximately 6 x 10 33 cm -2 . To proceed further, we Mellin transform equ. (^j), 
using also equs. @— (§). Defining 



h(r) 



J E k I(E, r)dE, 
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we obtain 



dl k 
dr 
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where 



Ck 
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= [jfj ( a c + a n - a n fk), 

h= [\l-y) k f(y)dy. 
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For k — 0, f = 1 and the total number of neutrinos is reduced because of the weak charged 
current, the weak neutral current providing only a redistribution of energy, fj. is a decreasing 
function of k, and for sufficiently large values of k, Ck corresponds to absorption created by 

Otot ( Otot = <?cc + One )• 

Let us define the column vector J n 

1 I ^ 
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lk+n/3 
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and the matrix 



where R n is the diagonal matrix 
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We obtain then 
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The formal solution of the above equation is 



J(r) = exp(Pr) J(r = 0). 



In terms of the first component 

1 - 2 1 

Jo = Jo + Ro Ji T + Ro R\ J2 t +■■■-] — : Ro Ri R2 ■ • ■ Rn- 
2! n\ 

with J k = J fc (r = 0). Finally in terms of the original moments we get 



hij) = h — Ck h+P T + f[C fc C k+ /3 4+2/3 r 2 H h 

+ (~^) n h. Ck Cfe +2/3 ' ' ' c k+(n-l)/3 h+n/3 T n + ■ ■ 



The above expression is the exact solution to equ. (|). 

The inversion of equ. in order to obtain I(E, r) is not an easy task. We prefer to 
study limiting behaviors of I(E, r) and then establish a consistent expansion. It is already 
stressed that fk tends to zero as k grows. At sufficiently large k, Ck becomes a constant 
independent of k 
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En 



(a c + a n ) =q 



(large k). 



(19) 



Substituting the above expression into equ. 



we get 



h(r) ~ J I(E)E k ex V [-a m {E)r] dE. (20) 

Therefore, at large energies (large k ) the neutrino intensity is well represented by the 
standard absorption formula 



labs 

(E,t) —1(E) exp [-a tot (E)T] 



(21) 



A systematic approach emerges, where the coefficients c& for k larger than some value are 
replaced by q (equ. (|19|)) and for the others we use the exact expression (equ. (0)). The first 
term in our expansion is obtained by keeping in equ. (|i~8D and replacing all Ck+ n p(n > 1) 
by q. We find 

-qr) a r 



n=0 



a c + a n 



n=0 
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(22) 



Returning to the energy variable we obtain 



h(E,r) = I abs (E,r) 



/-/, 



a c + a n ' 

Under the assumption of a power law initial neutrino spectrum 

1(E) = I(E ) ^ 

equ. (0) provides 



h(E,r)= 1(E) exp [-a tat (E)r) + 



a 6* 



(23) 



(24) 



1-exp -(^) a tot (E)r 



(25) 



Successive terms in our expansion can be obtained by using the exact expression for more 
than one Ck (two, three, ■ • •). 

One could treat the second term in equ. (|6|) as a source function F. Defining 

dy 



and multiplying equ. © by exp[a to t(E)r] results in 



y 



(26) 



— = exp (a tot r) F, 
dr 



(27) 



where H = exp(<T to tr)7. Direct integration yields 

I S (E, t) = 1(E) exp (-a tot r) + exp (-a tot r) f F(E, t) exp (a m t) dt. (28) 
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Reasonable approximations may be obtained by inserting "appropriate" forms for / in 
equ. (^). Assuming the generic form 



I P (E, t) = 1(E) exp [-a p (E)r] 
with cr p (E) to be specified, equ. ( pB"D provides 



(29) 



I s (E,t)= 1(E) exp [-cr tot r] + 







exp[— atot {E)t} -exp 


'-(^fME)r 




ap(E) (1 vV 1 





(30) 



Adopting Op = (i.e. I p = 1(E) ) we obtain an expression similar to equ. (|25|) . Both 
expressions are valid for short r and they overestimate significantly the actual result at 
large r. The other obvious choice is cr p (E) = cr tot (E), i.e. we insert in equ. (|26D the zeroth- 
order answer for I(E, r). We proceeded to a numerical solution of equ. (||) and in fig. 1 we 
show the shadowing factor defined by 

I(E,t) 



S(E, 



T 



1(E) 



(31) 



as a function of the energy for the maximum value of r. For the initial neutrino spectrum 
we used the form of equ. (p4|) with 7 = 2. Existing calculations of the unresolved AGN 



neutrino flux |T0|-[T2| give differential spectra with spectral index 7 varying between and 2 



over the energy range. The absorption factor suggested by equ. (|2T| ) is shown also in fig. 1 
(long dashed curve). It remains always below the numerical shadowing factor (solid line) 
and the two curves approach each other at large energies, as expected. The dash-dotted 
curve represents the shadowing factor implied by equ. (^) with u p (E) = <r tot (E). It is 
a significant improvement over the simple absorption formula (equ. (^l|)), although it lags 
behind the numerical result. 

Searching a more "appropriate" form for the neutrino intensity, we notice that for the 
specific form of f(y) (equ. (|5])), we have 

A=i -^i)( i+ ^4 + -+D' (32) 

Then for moderate values of k (k > 1, but not very large) c k remains relatively constant at 

( 1 \^ 

c k ~ I — (a c + < y > a n ) . (33) 



,E 

With the above c k , the neutrino intensity is provided by an absorption formula, equ 
with 

a p {E) = a cc (E)+ <y> a nc (E). (34) 



The shadowing factor given by equs. (|30| ) and ([34]) is shown in fig. 1 (short dashed curve) 
and it is in excellent agreement with the numerical simulation. 

Another approximate solution to the neutrino transport equation has been suggested in 
ref [0. Defining the effective cross section 

a eff (E) = cr cc (E) + a nc {E) - a nc {E) [' dy(l - yV^- l f{y), (35) 

Jo 

the neutrino "regeneration" effect is (incorrectly) exponentiated through the formula 

Ir(E, t) = 1(E) exp [-a eff (E)r] . (36) 

The ratio of the suggested shadowing factor by equ. (^) to the numerical one is shown in the 
inset of fig. 1 (dashed line). Clearly equ. (|36D overestimates the neutral current contribution. 
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In the same inset we show the corresponding ratio for our formula, equs. (|30| ) and ([34]) (solid 
line) . 

The neutral current removes neutrinos from the large energy part of the spectrum and 
populates the lower energy part of the spectrum. It is evident that the shadowing factor 
depends upon the initial spectrum, since a steeply falling spectrum does not provide consid- 
erable regeneration, while a spectrum extending to high energies (i.e. flat spectrum) induces 
large regeneration. Calculations of the AGN neutrino flux [IIJ suggest a flat energy distri- 



bution (7 = 0) up to energies E max = 10 3 TeV. Fig. 2 presents the shadowing factor for 
7 = with the other assumptions unchanged as in fig. 1. The disagreement of the shadowing 
factor implied by equ. (|36|) to the numerical calculation (dashed line in the inset) increases. 
Fig. 2 shows that, contrary to simple absorption, the shadowing factor at low energies can 
become greater than 1.0. This is an effect we could anticipate. Neutral current shifts the 
neutrino energy to lower values. If the new energy is low enough, such that atot{E)r <C 1, 
then the neutrino traverses the Earth unabsorbed. Therefore, at such nominal energy, apart 
from the initial neutrinos, we will have the neutrinos displaced in energy by the neutral 
current. In our calculations (fig. 2) the shadowing factor at 1 TeV is 1.36. 

The use of ultrahigh energy neutrinos in order to image the Earth's internal structure 



14Hl9| has been considered for some time. It is highly interesting to obtain further infor- 
mation on the density distribution of the Earth, independently of the seismological deter- 
minations. In these investigations the modified neutrino flux is given by simple absorption, 
ignoring the neutrino regeneration by the neutral current. We provided a simple, analytical 
and accurate description, equs. (^) and fl34|), for the neutrino propagation inside the Earth, 
for all the values of energy E and depth r. Our formalism will be essential for all reliable 
evaluations of the event rates in the neutrino telescopes under construction In the 

present work we used simple parametrizations for the neutrino cross sections, which may not 
be very accurate. Our main concern here is to establish a correct procedure, rather than to 
provide detailed numbers. At very high energy, the cross sections sample the nucleon parton 
densities at small x, a kinematical region where new physical phenomena (BFKL physics 
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pO| ) are operative. The shadowing factor is very sensitive to the actual magnitude of the 
cross section and therefore the neutrino passage through the Earth will provide information 
about the structure of the nucleon at small x. It seems that the detection of ultrahigh energy 
cosmic neutrinos will be important not only for astrophysics, but also for particle physics 
and geophysics. 

Acknowledgements. 

One of us (A.N.) would like to acknowledge useful discussions with P. Gondolo, V. 
Stenger, J. Learned, L. Resvanis, I. Sarcevic, F. Stecker, A. Tarantola, D. Kazanas and S. 
Ichtiaroglou. This research was supported in part by the Fulbright Foundation and the EU 
program "Human Capital and Mobility" . 



10 



FIGURES 



1 


■ i i i | 


1 


1 1 1 1 1 1 1 1 1 1 1 


_ 


0.8 




1 .4 
1 .2 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


- 






Ratio 

o o 

0) 00 - 1 






Shadowing Factor 
p " p 


It 
"ll 

u 

\ 1 







0.2 


\\ 

\V 

x V 

\ \- 




100 200 300 400 500 
Energy (TeV) 















1 


i 





100 200 300 400 500 

Energy (TeV) 

FIG. 1. The shadowing factor for neutrinos traveling along the Earth's diameter, as a function 
of the energy for spectral index 7 = 2.0. The long-dashed curve is the absorption factor given by 
equ.(21), the dot-dashed curve is the shadowing factor implied by equ.(30) with a p = atot, the 
short-dashed curve is the shadowing factor implied by equs.(30) and (34), while the solid line 
represents the numerical evaluation of the shadowing factor. The inset shows the ratio of the 
proposed shadowing factor to the numerical one. The dashed-curve is derived from equ.(36) and 
the solid curve is derived from equs.(30) and (34). 
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FIG. 2. The shadowing factor for neutrinos traveling along the Earth's diameter, as a function 
of the energy, for spectral index 7 = and E < 1100 TeV. Different curves correspond to the same 
conditions as in fig. 1. 
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